Detecting binary neutron star systems with spin in advanced gravitational-wave detectors. 
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The detection of gravitational waves from binary neutron stars is a major goal of the gravitational-wave 
observatories Advanced LIGO and Advanced Virgo. Previous searches for binary neutron stars with LIGO 
and Virgo neglected the component stars' angular momentum (spin). We demonstrate that neglecting spin in 
matched-filter searches causes advanced detectors to lose more than 3% of the possible signal-to-noise ratio for 
59% (6%) of sources, assuming that neutron star dimensionless spins, cJ/GM 2 , are uniformly distributed with 
magnitudes between and 0.4 (0.05) and that the neutron stars have isotropically distributed spin orientations. 
We present a new method of constructing filter banks for advanced-detector searches, which can create template 
banks of signals with non-zero spins that are (anti-)aligned with the orbital angular momentum. We show that 
this search loses more than 3% of the maximium signal-to-noise for only 9% (0.2%) of BNS sources with 
dimensionless spins between and 0.4 (0.05) and isotropic spin orientations. Use of this template bank will 
prevent selection bias in gravitational-wave searches and allow a more accurate exploration of the distribution 
of spins in binary neutron stars. 

PACS numbers: 04.80.Nn, 04.25.Nx, 04.30.Db, 



I. INTRODUCTION 

The second-generation gravitational wave detectors Ad- 
vanced LIGO (aLIGO) and Advanced Virgo (AdV) ffl E) 
are expected to begin observations in 2015, and to reach 
full sensitivity by 2018-19. These detectors will observe a 
volume of the universe more than a thousand times greater 
than first-generation detectors and establish the new field of 
gravitational-wave astronomy. Estimated detection rates for 
aLIGO and AdV suggest that binary neutron stars (BNS) will 
be the most numerous source detected, with plausible rates 
of ~ 40/yr [3|. Gravitational wave observations of BNS sys- 
tems will allow measurement of the properties of neutron stars 
and allow us to explore the processes of stellar evolution. 

The gravitational waves that advanced detectors will ob- 
serve from inspiralling BNS systems are well described by 
post-Newtonian theory Q. As the neutron stars orbit each 
other, they lose energy to gravitational waves causing them to 
spiral together and eventually merge. If the angular momen- 
tum (spin) of the component neutron stars is zero, the gravita- 
tional waveform emitted depends at leading order on the chirp 
mass of the binary M = (mirr^) 3 ^ 5 / (mi + m^ 1 ^ E), 
where mi, are the component masses of the two neutron 
stars, and at higher order on the symmetric mass ratio r\ — 
va\m,il(m\ + m-i) 2 Il6l- [TTI . If the neutron stars are rotating, 
coupling between the neutron stars' spin Si .2 and the orbital 
angular momentum L of the binary will affect the dynamics 
of BNS mergers lfT2HT31 . We measure the neutron stars' spin 
using the dimensionless parameter X1.2 = Si^/m\ 2 . 

The maximum spin value for a wide class of neutron star 
equations of state is \ = \x\ ~ 0.7 fl6l . However, the spins 
of neutron stars in BNS systems is likely to be smaller than 
this limit. The spin period at the birth of a neutron star is 
thought to be in the range 10-140 ms 02] QjD. During the 
evolution of the binary, accretion may increase the spin of one 
of the stars [19|, however neutron stars are unlikely to have 
periods less than 1 ms [20], corresponding to a dimensionless 



spin of x ~ 0.4. The period of the fastest known pulsar in a 
double neutron star system, J0737-3039A, is 22.70 ms HH. 
corresponding to a spin of only \ ~ 0.05. In this paper, we 
therefore consider two populations of neutron star binaries: 
the first has spins uniformly distributed from \ — to 0.4, 
the second, a sub-set of this, has spins between and 0.05. 
This extended spin distribution allows for the possibility of 
serendipitous discovery of BNS systems in globular clusters, 
where the evolutionary paths may be different than that in field 
binaries [22|. Since supernova kicks may cause the direction 
of the neutron star's angular momentum to be misaligned with 
the orbital angular momentum of the binary |23], or the bina- 
ries may be formed by direct capture, we consider a popula- 
tion of binaries with an isotropic spin distribution. 

Searches for binary neutron star systems in gravitational- 
wave detectors use template -based searches [24] . Data from 
the detector is correlated against a bank of known template 
waveforms, which cover the space of parameters searched 
over 0251 . The template bank is constructed so that it covers 
the parameter space of interest so that any signal in this region 
will lose no more than 3% of the signal-to-noise ratio obtained 
by an exactly matching template. Alternative search methods 
have been proposed [26 27], however these still require the 
construction of a template bank to perform the search. The 
effect of spin-orbit and spin-spin interactions were neglected 
in previous BNS searches [28 1, as they do not have a signif- 
icant effect on the ~ 1600 gravitational wave cycles in the 
40-2000 Hz sensitive band of first-generation detectors l29l . 
However, aLIGO and AdV will be sensitive to gravitational- 
wave frequencies between 10-2000Hz, increasing the number 
of cycles in band by an order of magnitude. Initial studies 
have demonstrated that over this band, the small secular ef- 
fects produced by spin-orbit and spin-spin coupling will have 
a significant effect on the detectability of BNS systems with 
non-trivial component spins ll30l . However, the current ge- 
ometric method for placing BNS templates OTl does not in- 
corporate spin. While numerical (stochastic) methods could 
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be used to include spin, these require substantially more tem- 
plates than a comparable geometric approach [32]. 

We present a new geometric algorithm for placing tem- 
plates for BNS systems with spin, which has a significantly 
higher sensitivity than previous searches. To demonstrate this, 
we first perform a systematic evaluation of the ability of a 
search that neglects spin to detect gravitational waves for BNS 
in aLIGO and AdV. We show that this search will lose more 
than 3% of the matched filter signal-to-noise ratio for 59% 
(6%) of signals if it is used to search for BNS systems with 
spins uniformly distributed between < \i,2 < 0.4(0.05); 
this is unsatisfactory over a large region of the signal param- 
eter space. We show that by considering BNS systems where 
the spin of the neutron stars are aligned with the orbital an- 
gular momentum (i.e. the binary is not precessing), we can 
create a two-dimensional template bank that is efficient at de- 
tecting spin-aligned BNS signals. Finally we demonstrate that 
this bank is sufficient to detect signals from generic spinning, 
precessing binaries in aLIGO and AdV. The spin-aligned bank 
loses more than 3% of the signal-to-noise ratio for only 9% 
(0.2%) of signals, sufficient to construct a sensitive and unbi- 
ased search for BNS systems in aLIGO and AdV. 



H. BNS SEARCH SENSITIVITY 

We quantify the performance of templated matched-filter 
searches by the fitting factor (FF) of the search [33 1. The 
fitting factor is the fraction of the signal-to-noise ratio that 
would be recovered when matching a given signal with the 
best matching waveform in the template bank. We first define 
the overlap between two templates hi and h% as 

which is defined in terms of the noise-weighted inner prod- 
uct ED 

WW = 4H«r». (2, 

J0 dn\J ) 

This overlap is the fraction of signal power that would be re- 
covered by searching for the signal hi using a matched filter 
constructed from hz- Maximizing the overlap over the time of 
arrival and waveform phase yields the match 

M(hi,h 2 ) =max.(h 1 \h 2 ((j) c ,t c )). (3) 

The mismatch, 1 — AL is the fraction of the optimal signal- 
to-noise ratio that is lost when searching for a signal hi with 
a template waveform h^- 

When searching for BNSs, we do not know the exact phys- 
ical parameters of the system. We assume that the masses 
of the neutron stars lie between 1 and 3 ill© and construct a 
bank of waveform templates to span this region of the 
mass parameter space. To measure the sensitivity of this bank 



to a gravitational waveform h s with unknown parameters, we 
compute the fitting factor 

¥¥{h s ) = max M(h s ,h), (4) 

he{h b } 

where we have maximized the match over all the templates 
in the bank. In searches for gravitational waves using LIGO 
and Virgo, the bank is constructed such that the fitting fac- 
tor for any signal in the target parameter space will never be 
less than 0.97. At least one of the templates in the bank must 
have a maximized overlap of 0.97 (or more) with the signal. 
This value is chosen to correspond to an event rate loss of no 
more than 10% of possible sources within the range of the de- 
tectors 1351. In this paper, we use a fitting factor of 0.97 to 
construct search template banks. 

We now test whether a bank of templates that does not 
model the effect of spin is sufficient to detect generic, spin- 
ning BNS sources in aLIGO and AdV. We create a bank of 
non-spinning templates that would recover any non-spinning 
BNS system with a fitting factor greater than 0.97. This 
bank is constructed using TaylorF2 waveforms, which are 
constructed using the stationary phase approximation to the 
gravitational-wave phasing accurate to 3.5 post-Newtonian 
(PN) order |4][36). To create a bank of these waveforms we 
use the hexagonal -placement method defined in [37], which 
was used in the majority of previous searches in LIGO and 
Virgo [38 -40 1 . This template bank is placed using the metric 
given in [25 1, which is valid, by construction, for templates at 
2PN order. Our signal waveforms are constructed using the 
SpinTaylorT4 waveform [41 1, a time-domain waveform accu- 
rate to 3.5PN order in the orbital phase which includes the 
leading order spin-orbit, spin-spin, and precessional modula- 
tion effects. We first confirm that although the bank is con- 
structed at 2PN order, it yields fitting factors greater than 0.97 
for both the TaylorF2 and SpinTaylorT4 non-spinning wave- 
forms at 3.5PN order. To simulate a population of spinning 
BNS sources, we generate 100,000 signals with component 
masses uniformly distributed between 1 and 3 Mq and di- 
mensionless spin magnitudes uniformly distributed between 
and 0.4. The orientation of the spin, the orientation of the 
orbital angular momentum, and the sky location are isotrop- 
ically distributed. To model the sensitivity of a second gen- 
eration gravitational wave interferometer, we use the aLIGO 
zero-detuned, high-power sensitivity curve ll42l . For our sim- 
ulations, we use a lower frequency cutoff of 15Hz. 

We note that for non-precessing systems the fitting factor is 
independent of the detector alignment and location; however 
this statement is not true for precessing systems. For such sys- 
tems, however, the distribution of fitting factors over a popula- 
tion of sources will be independent of the detector alignment 
and location. Therefore, for this study we calculate the fitting- 
factor for a single detector with an arbitrary location and po- 
sition. 

In Fig.[T]we show the distribution of fitting factors obtained 
when searching for our population of BNS sources with the 
non-spinning template bank. We see that 59% of signals were 
recovered with a fitting factor less than 0.97. If the maxi- 
mum spin magnitude is restricted to 0.05, we find that 6% of 
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Fitting Factor 

FIG. 1 . The distribution of fitting factors obtained by searching for 
the precessing BNS systems described in section [n] with component 
spins up to 0.4 (blue solid line), 0.2 (green dashed line), and 0.05 (red 
dotted line) using the non-spinning BNS template bank described in 
section [n] and the advanced LIGO, zero-detuned, high-power PSD 
with a 15Hz lower frequency cutoff. 



signals are recovered with a FF less than 0.97. If BNS sys- 
tems do exist with spin magnitudes up to 0.4, a template bank 
that captures the effects of spin will be required to maximize 
the number of BNS detections. Detection efficiency will be 
greatly reduced by using a template bank that only contains 
waveforms with no spin effects. Even under the assumption 
that component spins in BNS systems will be no greater than 
0.05, detection efficiency will be decreased if the effect of spin 
on the signal waveform is ignored. 



III. A TEMPLATE PLACEMENT ALGORITHM FOR 
ALIGNED-SPIN BNS TEMPLATES 

As we have demonstrated in the previous section, there 
is a substantial region of the BNS parameter space where a 
significant loss in signal-to-noise ratio would be encountered 



when searching for astrophysically plausible, spinning BNS 
systems with non-spinning templates. It has been suggested 
that using BNS templates where the spins of the system are 
aligned with the orbital angular momentum is sufficient for 
detecting generic BNS systems with second-generation detec- 
tors fTO | using TaylorF2 templates that incorporate the leading 
order spin-orbit and spin-spin corrections ||43l . 

In this section we use these spin-aligned waveforms to con- 
struct a template bank that attempts to cover the full space of 
astrophysically plausible BNS spin configurations. This tem- 
plate bank should contain as few templates as possible, while 
still being able to detect any BNS system that might be ob- 
served with aLIGO and AdV. To achieve this, it is important to 
assess the "effective dimension" of the space, which is defined 
as the number of orthogonal directions over which template 
waveforms need to placed in order to cover the full physically 
possible parameter range. While previous attempts to cover 
the BNS parameter space have made use of 3 or 4 dimen- 
sional template banks [44|, we demonstrate that the effective 
dimension of this parameter space is only two dimensional. 
For BNS systems in aLIGO and AdV the extent of the physi- 
cal parameter space in the remaining directions is smaller than 
the coverage radius of a template and can be neglected. 

As the effective dimension of the space is two-dimensional, 
a hexagonal placement algorithm, similar to that used in pre- 
vious searches of LIGO and Virgo data, could be employed to 
cover the space. This allows our new method to be incorpo- 
rated into existing search pipelines in a straightforward way. 

Since BNS systems coalesce at ~ 1500 Hz, significantly 
higher than the most sensitive band of the detectors, the wave- 
form will be dominated by the inspiral part of the signal B31 . 
The effect of component spin on BNS inspiral waveforms has 
been well explored in the literature lfT2lfl4l l46l ). For spin- 
aligned (i.e. non-precessing) waveforms, the dominant effects 
of component spin are spin-orbit coupling, which enters the 
waveform phasing at 1.5PN order, and spinl-spin2 coupling, 
which enters the waveform phasing at 2PN order. Other spin- 
related corrections to the PN phasing have been computed 
P7l l48l . however, in this work we restrict to only the two 
dominant terms. The methods described here are easily ex- 
tendable to include additional spin correction terms, however 
we do not expect them to significantly change our results. 



To construct a bank to search for generic BNS signals, we use TaylorF2 waveforms accurate to 3.5PN order in orbital phase 
and including the leading order spin-orbit and spin-spin terms given by Il43~ll45l 

h(f)=A(f;6 x y*^ (5) 

where 8 X describe the various orientation angles that only affect the amplitude and overall phase of the observed gravitational 
waveform E4ll . The phase W is given by 

* = 2Trf xt c ~(j> c + X^x-^ 3 + Xzx- 1 + X 3 x- 2 / 3 + A 4 ar 1/3 + A 5i log(x) + Agx 1 / 3 + X 6L \og{x)x 1 ^ + X 7 x 2 / 3 , (6) 



where / is the frequency, fo is a fiducial frequency, x — f / /q, t c is the coalescence time, cf) c is a constant phase offset. The PN 
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phasing terms are 
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where 75 is the Euler gamma constant, /3 (the dominant spin-orbit coupling term) and a (the dominant spin-spin coupling term) 
are given by 
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and L is the unit vector in the direction of the orbital angular momentum. Note that above we have omitted the A5 term, as it has 
no dependance on frequency and is therefore included in the constant phase offset, <j) c . 



Our goal is to construct a template bank containing the min- 
imum number of waveforms for which any plausible BNS sig- 
nal has a FF of 0.97 or higher. To place a template bank, we 
follow the method of Owen ll49l . We first construct a metric 
on the waveform parameter space that describes the mismatch 
between infinitesimally separated points, 

O(h(8),h{0 + 50)) =l-^2g ij (O)80 i 56 j , (17) 



with the metric given by, 
9iA e ) 
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times" To, T3 ll25l . In this coordinate system, ellipses are con- 
structed that describe fitting factors greater than 0.97 around 
a point and hexagonal placement is used to efficiently tile the 
space to achieve the desired minimal match l3D . 

To construct our new bank, we treat the six and two 
XiL components, given in Eq. Q, as eight independent pa- 
rameters, as in [50 1 . The range of possible physical values 
will trace out a four-dimensional manifold in the eight dimen- 
sional parameter space given by the X a , where a is an index 
that takes both i and iL values. We will demonstrate that this 
eight-dimensional parameter space allows us to construct a 
metric without intrinsic curvature. 



and where 9 describes the parameters of the signal, in this 
case the masses and the spins. 

We can use this metric to approximate the mismatch in the 
neighborhood of any point in the space by transforming to 
a coordinate system in which the metric is locally flat. To 
achieve this when placing the two-dimensional non-spinning 
bank, the masses mi,m2 are transformed into the "chirp 



As shown in [ 49 1 it is possible to solve equation ( 18 1, max- 
imizing over the phase, cj>c, to give the metric in terms of a 9 
dimensional space: 



(19) 
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In this expression tp a is given by 

a* 

ipo = wr = 27r/ a; 

<9A,; 



9* 

3A,r 



= a; (4 ~ 5)/3 log(x) 



(20) 
(21) 
(22) 



and J' is the moment functional of the noise PSD [ 43 , 49 1 
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and xu and ij, correspond to the lower and upper bounds 
of frequency in the integral. Unless stated otherwise we use 
h =x L f = 15Hz for the aLIGO PSD and choose 2000Hz 
for the upper frequency cutoff, fu — xu/q. While it is un- 
physical to use the same upper frequency cutoff for all sys- 
tems, especially as we are not including a merger in our wave- 
form model, it is necessary to make this assumption to ensure 
that our metric will be flat. For BNS systems this approxima- 
tion is fair to use as such systems will merge at frequencies 
that are outside the sensitive range of the advanced detectors 
and thus our calculation of signal power is not affected by as- 
suming that all BNS systems merge abruptly at 2000Hz. This 
approach was also used in |31 1 for computational efficiency. 
We do not go into detail here about the terms including log- 
arithms of the frequency, which first appear at 2.5PN order; 
these are discussed in detail in fBTl . 

Following [49] we can then maximize this expression over 
tc to give the metric in terms of the eight X a 



9a/3 = 
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(25) 



It is worth highlighting that the parameter space metric g a p, 
in the A Q coordinate system, has no dependence on the values 
of A Q . In other words, the parameter space is globally flat in 
this eight-dimensional parameter space. 

Although this eight-dimensional metric is globally flat, we 
have increased the dimensionality of the physical waveform 
space by a factor of two. However, we can transform this met- 
ric to a new coordinate system that will allow us to assess the 
effective dimensionality of the parameter space. We first rotate 
and rescale the metric to transform to a Cartesian coordinate 
system. We now use indicies i, j to number the remaining 
eight X a coordinates. As is a real, symmetric matrix we 
can use the eigenvalues and eigenvectors of the metric to ro- 
tate into an orthonormal coordinate system defined by 



(26) 



where Vy describes the eigenvectors of and its cor- 
responding eigenvalues. We use the convention that is 



the j th component of the i th eigenvector, and the eigenvec- 
tors are normalized by V T V = 1. In this coordinate system, 
the metric, g'^, will be the identity matrix. Next, we perform 
a rotation to align the axis of the parameter space with the 
principal components of the physically possible region of the 
space. The physically allowed ranges of the masses and spins 
cover only a limited region in the parameter space. The ex- 
tent of the physically relevant region of the space in a certain 
direction may be thin relative to the desired mismatch. By 
orienting the coordinate system along the principal directions 
we can easily identify any orthogonal directions in which the 
physical region is sufficiently thin that we do not need to place 
templates in those directions. This will allow us to assess 
the effective dimension of the parameter space, or, in other 
words, how many directions we need to consider when plac- 
ing a template bank. Transforming to a Cartesian coordinate 
system also helps with template placement, as it is trivial to 
place templates using the optimal A* n lattice ll52ll in a 2, 3 or 
4 dimension Cartesian coordinate system 

To perform the second rotation we make use of the fact that 
in a Cartesian coordinate system we are free to rotate the co- 
ordinates without changing the form of the metric. We would 
like to rotate the coordinates so that the greatest extent of the 
template bank lies along as few directions as possible. To ac- 
complish this we first draw many examples of physical param- 
eters of the masses and spins, and calculate the corresponding 
values of [ii for each of these points. We then do a Princi- 
pal Component Analysis on this dataset, which amounts to 
finding the eigenvectors of the covariance matrix from the set 
of fii. This produces a rotation into a new set of coordinates 
given by 



6 = £ fa 



(27) 



where CV, contains the eigenvectors of the covariance matrix 
using the same conventions as for Vij . The rotation of course 
leaves the metric Cartesian, but now the bank is oriented along 
the principal axes and it is much easier to visualize the shape 
of the boundaries and determine how to perform the template 
placement. 

We now use this method to construct a template bank 
where the spin of each component neutron star is restricted 
to 0.4. When this metric is constructed using the aLIGO, 
zero-detuned high-power noise curve with a lower frequency 
cut-off of 15Hz we show that, although many additional tem- 
plates are required to cover an aligned spinning parameter 
space when compared to the non spinning space, the effective 
dimension for these BNS systems is still two. 

We begin by attempting to visualize the space. We will refer 
to £i as the direction along which the parameter space has the 
biggest extent (the dominant direction) and £g as the direction 
with the smallest extent (the least-dominant direction). We 
draw a large set of points, with random values of masses and 
spins, and transform these points into the coordinate sys- 
tem. The position of these points is shown in Figure [2] where 
we plot the extent of £2,3,4 against £1. 

In Figure [2] and subsequent plots, we have scaled the £3 di- 
rection such that one unit corresponds to the coverage diam- 
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FIG. 2. The extent of the binary neutron star, \i < 0.4, parame- 
ter space in the £2, £3 and £4 directions, plotted against £1. The 
coordinates have been scaled such that one unit corresponds to the 
coverage diameter of a template at 0.97 mismatch. Generated using 
the zero-detuned, high-power advanced LIGO sensitivity curve with 
a 15Hz lower frequency cut off. 



0.9^J We remind that mismatch is proportional to distance 
squared and therefore two points separated by one unit would 
have a match of 0.88 

Immediately we notice that the extent along the £4 direc- 
tion is small compared to the diameter of a template. We can 
also see that the extent along the £3 direction is comparable 
to a template diameter, while the £1 and £ 2 directions have 
much larger extents and clearly need to be gridded over. The 
extent in the other 4 directions is smaller than £4 and can be 
completely ignored. This hierarchy of measurable parameters 
may be a generic feature according to the model of 1 53 1 . 

The plot of £1 against £3 in Figure [2]can be somewhat mis- 
leading as we have projected out the £2 direction. It is more 
informative to investigate the depth of £3 at fixed values of 
£1 and £ 2 an d translate this into the maximum mismatch that 
would be obtained if one were to assume that there is no width 
in the third direction. In Figure|3]we show the maximum mis- 
match between the central and extremal values of the possi- 
ble range of £3 (and £4) as a function of the two primary di- 
rections. This is calculated by binning the points mentioned 
above into bins in £1 and £2, where the bin width is equal to 
one template radius. We then determine the extremal values 
of £3 (and £4) for the points in each bin. From Figure[3]we can 
see that, while there are small areas of parameter space where 
up to a 1.6% loss in SNR would be incurred from assuming the 
£3 direction had no depth, most areas of the parameter space 
are very thin in the £3 direction. This figure also helps to 
reinforce the fact that the depth in the fourth direction is neg- 
ligible, as, even in the worst region of the space, no more than 
0.01% of SNR would be lost by assuming £ 4 had no depth. 
The depth of the £5 3 directions are even smaller than £4. 

In this coordinate system it is easy to explore how the size 
of the parameter space depends on the maximal spins of the 
component neutron stars. In Figure |4] we show the extent of 
the physical space for aligned spinning BNS systems, with 
maximum component spins of 0.4, 0.2 and 0.1, compared to 
that of non-spinning systems. Ignoring any issues related to 
the depth of the £ 3 direction, one can clearly see that to cover 
the aligned spin parameter space will require a great deal more 
templates than the non spinning parameter space. 

From these results we can see that a 2 dimensional template 
bank would be sufficient to cover the aligned spin parameter 
space for BNS systems in the advanced detector era. Specifi- 
cally, we would advocate placing a hexagonal lattice in the £1, 
£2 coordinates and setting the value of £ 3 8 to be the middle 
of the possible range of those parameters at the given position 
of £1, £2- For the regions of parameter space where the depth 
of £3 is not negligible, one could either ignore it, understand- 
ing that the resulting bank will not have a fitting factor of 0.97 
in this region. Alternatively, one could stack templates in the 
region where £3 is deepest to minimize this effect. 

For this work we chose to employ a hexagonal template 
bank in the £x> £2 coordinates, stacking the templates in the 
£3 direction, where necessary, to ensure that the maximum 



eter of a template at 0.97 mismatch. Equivalently, we have 
scaled the directions such that two points separated by 0.5 
units (one template radius) in any direction have a match of 



The unsealed distance between two points with a match of 0.97 would be 

(1 - 0.97) 05 = 0.17 
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FIG. 3. The mismatch between the edge and centre of the physically 
possible range of £3 (top) and £4 (bottom) values as a function of 
£1 and £2. The fi coordinates have been scaled such that one unit 
corresponds to the coverage diameter of a template at 0.97 mismatch. 
Plotted for a binary neutron star parameter space with spins restricted 
to 0.4 using the zero-detuned, high-power advanced LIGO sensitivity 
curve with a 15Hz lower frequency cut off. 



mismatch due to the depth of £3 is no more than 0.25%. For 
an aligned-spin template bank where the spin of each com- 
ponent is restricted to 0.4, using the advanced LIGO, zero- 
detuned high-power noise curve with a lower frequency cut- 
off of 15Hz, we find that approximately 520,000 templates are 
required. Roughly 100,000 of these templates were added by 
the stacking process. 

We can verify that the template bank algorithm is working 
correctly by repeating the simulation described in section |U] 
but evaluating the fitting factor between our bank of aligned- 
spin template waveforms and a set of signals that is restricted 
to having spins that are (anti-)aligned with the orbital angular 
momentum. The results of this simulation are shown in figure 
[5]and one can see that with our bank we do not observe fitting 
factors lower than 0.97 when searching for aligned spin BNS 
systems. 

In the previous paragraphs we have restricted attention 
to the aLIGO zero-detuned, high-power predicted sensitivity 
with a 15Hz lower frequency cut off. However, we should 
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FIG. 4. The size of the BNS parameter space as a function of the 
maximum spin. The darkest points indicate points with spin on both 
components constrained to 0.4, then, in order of increasing lightness, 
we show points constrained to a maximal spin of 0.2 and 0.1, finally 
the lightest points show points with no spin. The coordinates have 
been scaled such that one unit corresponds to the coverage diame- 
ter of a template at 0.97 mismatch. This plot was generated using 
the zero-detuned, high-power aLIGO sensitivity curve with a 15Hz 
lower frequency cut off 
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FIG. 5. The distribution of fitting factors obtained by searching 
for aligned spin, binary neutron star systems, with spin magnitudes 
restricted to 0.4 using the aligned-spin BNS template bank described 
in section [Til] and the aLIGO, zero-detuned, high-power PSD with a 
15Hz lower frequency cutoff. 



verify that the conclusions we have drawn are valid for AdV, 
whose PSD is different from that of aLIGO, as shown in Fig- 
ure^ Additionally we should also show that the choice to use 
a 15Hz cut off in the aLIGO PSD does not affect the conclu- 
sions made in this section. 

The process we described above is applicable for any PSD, 
and therefore we can use it directly to determine the £j di- 
rections for the AdV PSD, or the aLIGO PSD with a 10Hz 
lower frequency cutoff. In Figure [7] we plot £1 against £ 2 for 
both PSDs while the color shows the mismatch between the 



8 



10" 



io- 



io- 



> 



2 

o 



10" 



10- 



10 



Advanced LIGO 
Advanced Virgo 
Initial LIGO 
Initial Vinro 



10 2 10 3 

Frequency (Hz) 



FIG. 6. The amplitude spectral density for the aLIGO zero-detuned 
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centre and edges in the £3 direction. This plot can be directly 
compared to Figure [3] We notice that the size of the param- 
eter space for the AdV PSD is significantly smaller than for 
the aLIGO PSD in all 3 of the dominant directions. Therefore 
our conclusions for aLIGO are still valid for AdV. Using our 
method we find that we require approximately 120,000 tem- 
plates to cover the parameter space for AdV, in comparison to 
approximately 520,000 templates for aLIGO. 

By comparing the results when using the aLIGO PSD with 
a 10Hz and 15Hz lower cut off we observe that using a 10Hz 
lower frequency cut off will increase the number of necessary 
templates from — 520000 to ~ 860000. However the shape 
of the parameter space, and thus our final conclusions, are un- 
affected when using a 10Hz lower frequency cutoff. However, 
in this case we see larger mismatches due to the depth of £3 
and therefore the process of stacking templates is important 
when using a 10Hz lower cut off. However, even in this case, 
we do not feel that the depth is large enough everywhere in 
the space to justify using a fully 3-dimensional placement al- 
gorithm. 



IV. COMPARISON TO ALTERNATIVE PLACEMENT 
METHODS 

An alternative approach to template placement for aligned 
spin systems is to use templates with "unphysical" values of 
the symmetric mass ratio, 77. That is, to use non-spinning tem- 
plates, with the desired range of chirp mass but where the 
range of 77 values is extended to include both values of r\ that 
are much lower than the relevant parameter space and values 
of 77 that are much higher, including templates with 77 greater 
than the physically possible limit of 0.25. 

We can understand this unphysical r\ approach in terms of 
our £j coordinate system by noting that it is always possible to 
produce a template with any possible value of £1 and £ 2 that is 
within the BNS parameter space, by using non-spinning tem- 
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FIG. 7. The mismatch between the edge and centre of the third 
dominant direction as a function of the first and second dominant 
directions when using the Virgo noise curve (top) and when using 
the advanced LIGO noise curve with a 10Hz lower frequency cut 
off (bottom). The £i coordinates have been scaled such that one unit 
corresponds to the coverage diameter of a template at 0.97 mismatch. 
Plotted for a binary neutron star parameter space with spins restricted 
to 0.4. 



plates with unrestricted values of 77. By generating a set of 
templates in the £1, £2 directions, where we restrict the chirp 
mass to be that possible for BNS systems, but where 77 ranges 
from 0. 1 to 0.7 we are able to cover the full physically pos- 
sible space in £1, £ 2 - However, the disadvantage to using un- 
physical 77 templates is that the points will not take the correct 
values of £3. The colorbar on Figure [HJindicates the mismatch 
between unphysical 77 templates and aligned-spin templates as 
a function of £1 and £2- In making this plot we assume that 
£3 has no depth in the aligned spin case by taking the central 
value where £3 has a range of values. 

While unphysical 7/ templates will produce an increase in 
efficiency when compared with non-spinning templates, the 
method is not as efficient as the aligned spin geometrical 
placement we have described. In addition, both methods re- 
quire the same number of templates to cover the parameter 
space. Therefore, we would recommend using aligned spin 
templates placed using our metric algorithm as opposed to un- 
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In making this plot we assume that £3 has no depth in the aligned 
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This plot was generated with spins restricted to 0.4 using the zero- 
detuned, high-power advanced LIGO sensitivity curve with a 15Hz 
lower frequency cut off. 



physical rj templates. 

Finally, we wish to compare the performance of this ge- 
ometrical algorithm with the stochastic bank proposed in 
I32ll54l . The stochastic placement works by randomly plac- 
ing points within the parameter space and rejecting points that 
are too "close" to points already in the bank. This has the ad- 
vantage that it is valid for any parameter space metric, so we 
could use any of the metrics discussed above. However, it is 
more computationally efficient to use the Cartesian & or 
coordinate system rather than the non-Cartesian metric given 
above. 

The disadvantage to a stochastic bank, when compared to 
a geometrically placed bank, is that it will require more tem- 
plates to achieve the same level of coverage 11321 . For our 
parameter space, consisting of BNS signals with component 
spins up to 0.4 and using the advanced LIGO zero-detuned 
high-power design curve with a 15Hz lower frequency cut-off, 
we found that the stochastic placement produced a bank con- 
taining ~ 750000 templates, which is 44% more than with the 
geometrical placement. However, stochastic placement can 
still be used to place templates when no analytical metric is 
known, such as when the merger becomes important. In such 
regions of parameter space, the stochastic placement may still 
be the best algorithm to use to place a template bank [44 1. 



V. PERFORMANCE OF THE ALIGNED SPIN TEMPLATE 
BANK 

In this section we would like to investigate the improve- 
ment in the detection of generic BNS systems that results 
from using a template bank that includes the dominant, non- 
precessing, spin effects. To do this we use the aligned spin- 
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FIG. 9. The distribution of fitting factors obtained by searching for 
the precessing signals described in section [n] with component spins 
up to 0.4 (blue solid line), 0.2 (green dashed line), and 0.05 (red 
dotted line) using the aligned spin BNS template bank described in 
section [Hi] and the advanced LIGO, zero-detuned, high-power PSD 
with a 15Hz lower frequency cutoff. 



results of using a nonspinning bank as shown in section [H] 

Using our aligned spin template bank, we repeat the in- 
vestigation from section [II] We create a population of source 
BNS signals identical to those used in [TTJ and compute the 
fitting factor between these signals and the aligned spin tem- 
plate bank. The results of this are shown in FIG|9] To de- 
crease the computational cost of this test, we only calculated 
the overlaps between a signal and templates that were within 
a range of ±0.1M Q in chirp mass. This is reasonable because 
the overlap will decrease rapidly with small changes in chirp 
mass, therefore we expect templates with very different val- 
ues of chirp mass to have low overlaps with each other. We 
verified that this approach did not cause us to underestimate 
the fitting-factor of our banks. 

We can now compare the results obtained in this section, us- 
ing our aligned-spin template bank, with the results obtained 
in section |IIJ using a non-spinning template bank. One can 
clearly see an improvement in the distribution of fitting fac- 
tors when using the aligned spin template bank. The fraction 
of signals that fall below a fitting factor of 0.97, when the 
spin magnitudes are restricted to 0.4, falls from 59% to 9%. 
We also see an improvement for signals that have spin mag- 
nitudes restricted to 0.05, where the fraction of signals falling 
below a fitting factor of 0.97 drops from 6% to 0.2%. We 
can also compare the performance of the aligned-spin bank to 
that of the non-spinning bank as a function of the maximum 



ning bank that we detail in section III and compare this to the 



spin magnitude, as shown in Figure 10 From this Figure we 
can see that regardless of the maximum component spin, the 
aligned spin bank will greatly reduce the number of signals 
recovered with fitting factors less than 0.97. 

A small fraction of signals fall below a FF of 0.97, even 
when using the new aligned-spin template bank. We expect 
that these poor matches with the aligned template bank are 
due to precession. In general, precessional effects will not be 
important in BNS systems as the orbital angular momentum 
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FIG. 10. The fraction of the precessing signals described in section 
[n] recovered with a fitting factor less than 0.97 as a function of the 
maximum component spin. Shown for the non-spinning BNS tem- 
plate bank described in section [TT] (blue solid line), and the aligned 
spin BNS template bank described in section [ill] (red dotted line). 
The advanced LIGO, zero-detuned, high-power PSD with a 15Hz 
lower frequency cutoff was used when computing the fitting factors. 

is significantly larger than the component spins. In such cases 
there is only a small angle between the total and orbital an- 
gular momenta and precession has only a small effect on the 
waveform. 

However, there is a small region of parameter space where 
precessional effects will have an effect for BNS systems. Us- 
ing the model of Ref. Il55l . applied to the small precession 
angles in BNS systems, we can predict for which systems pre- 
cession will be most important. The orientation of a precess- 
ing binary must be defined using the total angular momentum 
rather than the orbital angular momentum as done with non- 
precessing binaries. The orientations with the worst matches 
should be those where the system is edge-on (angular momen- 
tum perpendicular to the viewing direction) and where the de- 
tector is nearly insensitive to the plus polarization and only 
sees the cross polarization (a binary overhead of the detector 
would have its angular momentum oriented 45° between the 
arms of the detector). We find that this is indeed the case; in 
fact, all cases with fitting factors less than 0.95 are close to 
this configuration. All of these cases also have biases in the 
recovered mass and spin parameters due to the secular effects 
of precession on the phasing of the waveform. 

VI. CONCLUSION 

In this work we have investigated the effects of neglect- 
ing spin when searching for binary neutron star systems in 



aLIGO and AdV. We have found that, if component spins in 
binary neutron star systems are as large as 0.4, then neutron 
star spin cannot be neglected, and there is a non-trivial loss 
in signal-to-noise ratio even if the maximum spin is restricted 
to be less than 0.05. We have developed a new algorithm for 
placing an aligned spin template bank in the BNS parame- 
ter space. We have shown that this bank works for aligned 
spin systems and have demonstrated that it does significantly 
better for generic, precessing BNS systems than the tradi- 
tional non-spinning bank. However, for the BNS aligned spin 
Xi < 0.4 parameter space the aligned spin bank requires ap- 
proximately five times as many templates as the non-spinning 
bank. This increased number of templates will increase the 
computational cost of the search and increase the number of 
background events, so needs to be balanced against the poten- 
tial gain in being able to cover a larger region of parameter 
space. A further advantage of our method is the ease with 
which it can be incorporated into existing or future search 
pipelines, which include the use of signal-based vetoes ll56l 
and coincidence algorithms lf57l . In future work we will in- 
vestigate how this template bank performs in data from the 
aLIGO and AdV detectors which includes non-Gaussian and 
non-stationary noise features. Finally we note that the method 
proposed in this work should be applicable wherever the Tay- 
lorF2 waveforms closely represent actual gravitational wave- 
forms. In a future work we will investigate how well this 
method performs in the binary black hole and neutron-star, 
black-hole regions of the parameter space. Wherever the Tay- 
lorF2 approximation begins to break down, a stochastic bank 
placement may still be the most viable option. 
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